Suppression of small scale dynamo action by an imposed magnetic field 
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Non-helical hydromagnetic turbulence with an externally imposed magnetic field is investigated 
using direct numerical simulations. It is shown that the imposed magnetic field lowers the spec- 
tral magnetic energy in the inertial range. This is explained by a suppression of the small scale 
dynamo. At large scales, however, the spectral magnetic energy increases with increasing imposed 
field strength for moderately strong fields, and decreases only slightly for even stronger fields. The 
presence of Alfven waves is explicitly confirmed by monitoring the evolution of magnetic field and 
velocity at one point. The frequency ui agrees with v\k\, where va is the Alfven speed and k\ is 
the smallest wavenumber in the box. 



I. INTRODUCTION 

Turbulent magnetic fields are seen in many astrophys- 
ical settings 0, 0, Q • Such magnetic fields usually result 
from the conversion of kinetic energy into magnetic en- 
ergy, i.e. from dynamo action. Numerical simulations 
show that a dynamo-generated magnetic field can be of 
appreciable strength even when there is no kinetic he- 
licity 0, [j| . Simulations have recently also shown that 
at scales smaller than about five times the energy carry- 
ing scale the magnetic energy spectrum seems to enter 
an inertial subrange where the magnetic spectral energy 
exceeds the kinetic spectral energy 0. This means that 
over any subvolume, whose scale is within the inertial 
range, there is always a larger scale component of the 
field with significant strength. This raises the question 
whether one can model the small scale properties of such 
turbulence simply by imposing a magnetic field. 

A lot of work has already been devoted to studying 
hydromagnetic turbulence in the presence of an exter- 
nal field ]3, H, El • Nevertheless, the super-equipartition 
magnetic energy seen in simulations without imposed 
field has never been seen in simulations with imposed 
field. An exception is when the magnetic Prandtl num- 
ber is large|lOj. However, the super-equipartition is then 
seen between the viscous and the resistive cutoff - not 
in the inertial range. It is one of our goals to elucidate 
this puzzle. Likewise, although dynamos with helicity 
can produce substantial super-equipartition on the scale 
of the system, they too are not able to produce super- 
equipartition in the inertial range In that sense the 
difference between dynamos with and without imposed 
field is similar to the difference between helical and non- 
helical dynamos. 

The views on the effects of external fields are divided. 
A common scenario that applies when the conditions for 



dynamo action are not met (e.g. if the magnetic Reynolds 
number is too small), is one where a local magnetic field 
can be enhanced simply by winding up an external mag- 
netic field. Possible candidates where this may be the 
case are Io and Ganymede, in which convection interacts 
with the field of Jupiter leading to local field enhance- 
ment (l2l Il3l Hil ] . A similar possibility may also apply 
to the solar convection zone where the large scale field of 
the 11 year solar cycle is primarily located at the bottom 
of the convection zone [l5j , but the overlying convection 
zone may shred the field to produce small scale field . 
Another possibility that has been discussed more recently 
is that the small scale field at the solar surface could be 
generated locally by a small scale dynamo operating near 
the surface |l7j . 

In hydromagnetic turbulence theory magnetic and ki- 
netic energy are assumed to cascade from large to small 
scales, similar to the hydrodynamic case, although re- 
cent work has established a strong intrinsic anisotropy 
[Tc| , which has no counterpart in the hydrodynamic case. 
However, this theory does not address the possibility of 
dynamo action. It remains therefore an open question as 
to what is the nature of the interaction resulting from 
imposed and dynamo-generated magnetic fields. In par- 
ticular, we shall present evidence that the imposed mag- 
netic field does not enhance dynamo action. Instead, the 
external field does actually suppress dynamo action, al- 
beit in a subtle way because the rms turbulent velocity 
is generally not decreased by a modestly strong magnetic 
field. We show that the suppression can be associated 
with the work term resulting from the Lorentz force due 
to the imposed field. It turns out that this term changes 
sign above a certain field strength such that a certain 
fraction of magnetic energy flows backwards to enhance 
the kinetic energy instead. 



II. EQUATIONS 
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We adopt an isothermal equation of state with con- 
stant (isothermal) sound speed c s , so the pressure p is 
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TABLE I: Summary of the different runs with forcing at 
fcf = f .5. All runs have magnetic Prandtl number unity. The 
field strengths 0.06, 0.3, and 3.0 correspond roughly to 0.5, 
2.0, and 20 times B cq = ^/Jtopou TlnB . 



Run 


Res. 


V 


= V 


ReM 


^rms 


^rms 


B 


B4 


128 3 


4 x 


1Q- 4 


280 


0.076 


0.17 


0.3 


CI 


256 3 


2 x 


w- 4 


400 


0.062 


0.12 





C2 


256 3 


2 x 


io- 4 


400 


0.070 


0.12 


0.01 


C3 


256 3 


2 x 


w- 4 


370 


0.094 


0.12 


0.06 


C4 


256 3 


2 x 


io- 4 


500 


0.088 


0.19 


0.3 


C5 


256 3 


2 x 


w- 4 


500 


0.075 


0.15 


3 


D4 


512 3 


1 X 


w- 4 


930 


0.089 


0.14 


0.3 


Ef 


1024 3 


8 x 


10~ 5 


1000 


0.075 


0.12 






related to the density p by p = pc^. The equation of 
motion is written in the form 



Dm 
~Dt 



-c 2 s Vlnp- 



J x B 



(1) 



where D/Di = d/dt + u ■ V is the advective derivative, 
J = V x B / fiQ is the current density, po is the vacuum 
permeability, 

F visc = v (V 2 « + |VV • u + 2S • Vlnp) , (2) 

is the viscous force, v = const is the kinematic viscosity, 



duj 
Ox, 



9u j 



-5 l3 V ■ u 



(3) 



is the traceless rate of strain tensor, and / is a random 
forcing function that consists of non-helical plane waves; 
see Refs jy, [19| for details. The continuity equation is 
written in terms of the logarithmic density, 



Dlnp 
Bt 



= -V • u, 



(4) 



and the induction equation is solved in terms of the mag- 
netic vector potential A, 



dA 
~dt 



ux B + r/V 2 A, 



(5) 



where r\ — const is the magnetic diffusivity, and B = 
B + b is the magnetic field consisting of the imposed 
uniform (fc = 0) field, B , and the deviations from the 
imposed field b = V x A. This split is necessary because 
the vector potential corresponding to Bo cannot be pe- 
riodic, while both B and A can well be assumed to be 
periodic. 

In the simulations summarized in Table |I] we have used 
the same method as described in Ref. [l9|. Kinetic and 
magnetic Reynolds numbers are defined as 
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FIG. 1: Sketch of the energy budget showing the kinetic 
and magnetic energy reservoirs together with the flow of en- 
ergy for Run D4. The numbers in parentheses correspond to 
Run El without an imposed field. The total dissipation rate 
is denoted by e = ek + cm, which is the sum of kinetic and 
magnetic energy dissipation rates. The + sign on Wlo and 
the direction of the corresponding arrow emphasize that, at 
least for sufficiently strong imposed fields, energy flows from 
the magnetic to the kinetic energy reservoir. 



respectively. Here, k[ is the average forcing wavenumber 
and PrM = v jr\= ReM /Re is the magnetic Prandtl num- 
ber. In all cases studied below we assume PrM = 1- We 
study cases where kf is either 1.5 or 5. 

The Pencil Code [2(j is used for all our simulations. 
The resolution is varied between 128 3 and 1024 3 mesh- 
points. Although we solve the compressible equations, 
the sound speed is large compared with the turbulent ve- 
locities. We find that the energies of solenoidal and po- 
tential components of the flow have a ratio E pot /E so i « 
10 _4 -10 -2 for most scales; only towards the Nyquist fre- 
quency the ratio increases to about 0.1. Thus, our results 
should be close to the incompressible limit. 

We use non-dimensional quantities by measuring 
length in units of 1/fci (where ki = 2ir/L is the small- 
est wavenumber in the box of size L; in the present case 
L = 2tt), speed in units of the isothermal sound speed 
c s , density in units of the initial value po, and magnetic 
field in units of (poPoc 2 ) 1 / 2 ■ 



III. ENERGY BALANCE 

A sketch of the overall energy budget is given in Fig. 
where we show the magnetic and kinetic energy reser- 
voirs together with arrows indicating the flow of energy. 
The arrow pointing into the kinetic energy reservoir is the 
energy flux e entering the simulation through the exter- 
nal forcing, while the arrows pointing to the right from 
the kinetic and magnetic energy reservoirs denote vis- 
cous and Joule dissipation, i.e. ck and cm, respectively. 
On the average and in the statistically steady state we 
expect e = ck + cm- 

The two arrows between the kinetic and magnetic en- 
ergy reservoirs correspond to the contributions to the 
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work done against the Lorentz force. In general this work 
term can be written as (u ■ (j x B)), where we have used 
j = J to emphasize that the current density has van- 
ishing volume average. However, since B = Bq + b, 
where b — V x A is the departure from the imposed field 
(also with vanishing volume average), we can divide the 
work term into a contribution from the fluctuating field, 
(u-(jxb)), and one from the imposed field, (it- (j x Bq)) . 
The latter can also be written as Bq ■ (ux j), which em- 
phasizes the fact that this term is quadratic in the fluctu- 
ations and can hence only transfer energy between kinetic 
and magnetic energy reservoirs at the same wavenumber. 
We can thus write 



dE 



M 



df 



-W L - W L o - £m, 



(7) 



where Em = (b 2 ) / /iq is the magnetic energy (per unit 
volume) of the induced field without the imposed field, 
Wl = (u ■ (j x b)) is the work done by the fluctuating 
fields, Wlo = Bo-(uxj) is the work done against winding 
up the imposed mean field, and 6m = WoC?' 2 ) is the loss 
from Joule heating. In a closed or periodic system such as 
the one considered here, there are no surface terms, which 
is why there is no term associated with the Poynting flux 
in Eq. 0. 

The numbers on the arrows give the energy fluxes 
for a simulation with a moderately strong imposed field 
(Run D4). The numbers in parentheses are the corre- 
sponding values for a simulation without imposed fields. 
By comparing the two we see that with an imposed field 
the content of the magnetic energy reservoir is slightly 
increased. Nevertheless, magnetic dissipation has de- 
creased and kinetic dissipation has increased. This sug- 
gests that an imposed magnetic field quenches the dy- 
namo. 

Naively one might have expected that the — Wlo term 
always "helps" the dynamo and that it therefore al- 
ways transfers energy from kinetic to magnetic energy 
by winding up the imposed field. This is however not 



FIG. 3: Contribution to the spectral energy transfer between 
kinetic and magnetic energies due to the imposed field. Power 
spectrum of Bo • {uk x jj.) normalized by Ejf for run C4 (solid 
line) and run C3 (dashed line). We clearly see that at the box 
scale there is transport of energy from kinetic to magnetic 
field, while at all other scales the transport is in the opposite 
direction, i.e. there is a suppression of the magnetic field. 
It is also clear that the suppression is much stronger at the 
smallest scales. 



the case. In Fig. [21 we show Wl and Wlo > normalized 
by £o — fcfPo u ?msi as functions of imposed field strength. 
(In those units the total energy input to the system is 
e sa 0.07eo) The negative contribution from —Wlo f° r 
large field strengths is actually the main reason that sim- 
ulations with strong imposed fields have less magnetic en- 
ergy; see Table HI Since the Who term is local in A:-space 
it also explains the general increase in kinetic energy at 
all scales. 

To quantify the above statement, we discuss now the 
spectral energy transfer function, Wlo(&) = B -(ukXjk), 
where Uk and jk are Fourier filtered velocity and current 
density. In Fig. ||we plot the ratio of W^oik) to the 
magnetic energy spectrum Eyi{k), divided by the eddy 
turnover time r = (k{U lmB )~ , using data from Runs C3 
and C4. Here, the spectra are normalized such that 
JW w (k)dk = W w and fE M {k)dk = E M . It turns 
out that, first, Wi,o(k) has a positive contribution to the 
magnetic energy at small wavenumbers. This explains 
the increase in magnetic energy at large scales (small 
wavenumbers). Second, at moderate and large wavenum- 
bers, Who(k) is positive, which explains the suppression 
of the magnetic energy. 



IV. SPECTRAL ENERGY CHANGES 

Next we investigate the effect of varying the strength 
of the imposed field on the magnetic and kinetic energies 
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FIG. 4: Magnetic (top) and kinetic (bottom) spectral energy 
at wavenumbers 2 and 32 as a function of Bo. fef = 1.5. 



FIG. 5: Same as Fig. |I]but with forcing at fef = 5. 



at different wavenumbers; see Fig.Q] We see that around 
Bq « B cq the magnetic energy is somewhat enhanced at 
small wavenumbers (k = 2, corresponding to modestly 
large scales), but decreased at large wavenumbers (k = 
32, corresponding to small scales). At the same time the 
effect on the velocity field is weak, but there is generally 
a tendency for enhanced velocities, especially at large 
scales. 

When the forcing is at kf — 5, instead of at kf = 1.5, 
the trends are very similar as in Fig. 0J see Fig. |3J In 
particular, at large scales there is first an increase and 
then a decrease of the magnetic energy as the imposed 
field strength is increased, while for small scales the mag- 
netic energy decreases for all imposed field strengths. We 
therefore conclude that the suppression of the magnetic 
field is, at least qualitatively, independent of the forcing 
scale. 

The situation is different in the presence of helicity 
where it has been argued that it is particularly the large 
scale magnetic field at k = ki that is affected by the in- 
troduction of an imposed magnetic field 21]. This can 
be interpreted as a suppression of the a effect. Repeating 
the simulations of Ref. l21l . we were able to confirm their 
findings; see also Ref. |22j . We also find that kinetic and 
magnetic energy spectra fall almost on top of each other 
when Bq = 2B ea . This is just like in the case without 
imposed field [llj, except at k = k\ where there is an 
additional field component due to the a effect. Further- 



more, increasing the field to Bq = 20 B oq we do recover 
the same suppression of the dynamo as without helicity, 
i.e. the spectra look similar to those of Run C5. Thus, 
the suppression of dynamo activity by the imposed field 
is rather general and affects equally helical and nonhelical 
dynamos. 

It is generally believed that hydromagnetic turbulence 
can be described as an ensemble of Alfven waves. This is 
true both for the Goldreich-Sridhar [18| and Iroshnikov- 
Kraichnan |U|^] theories. This would then suggest that 
magnetic and kinetic energies should be comparable to 
one another at each scale. From Figs. El and we see that 
magnetic and kinetic energies are close to each other, but 
generally not equal. This is also seen in the simulations 
of Cho & Vishniac Only when the imposed field is 
approximately equal in strength to the rms field do we 
have approximate equipartition between magnetic and 
kinetic energies at small scales. 



V. SHAPE OF THE ENERGY SPECTRA 

As the resolution is increased, one begins to see indi- 
cations of the build-up of a short fc~ 5 / 3 inertial range of 
kinetic and magnetic energies at intermediate wavenum- 
bers; see Fig. |S1 The inertial range is as yet too short to 
be conclusive, and we therefore need larger simulations 
in order to be sure whether we have a real fc~ 5 / 3 slope or 
not. 
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FIG. 6: Magnetic and kinetic energy spectra for runs with FIG. 7: Magnetic and kinetic energy spectra for runs with 
different imposed field strengths (Runs CI and C3-C5). In all different imposed field strengths and forcing at k — 5. 
cases B aq = 0.12-0.15; see Table □ 



so 



From Fig. [S] we also see that in the range fci < k < 
10 the magnetic energy spectrum seems to follow a k~ l 
slope. For comparison, in the case without an imposed 
field the spectral magnetic energy was actually increasing 
with k and followed approximately a fc 1 / 3 slope 0] at 
small k. The k~ l spectrum for imposed fields can be 
motivated by dimensional arguments: assume that the 
magnetic energy spectrum is a function of the imposed 
field strength Bq and the wavenumber k, and that the 
spectrum is given by the ansatz Eu(k) — CB^k h , then, 
from dimensional arguments, one finds a = 2 and b = — 1, 



E M (k) = CBlk- 1 (8) 

where C is a dimensionless constant. Such a spectrum is 
expected if there is a mean field |25j . but it may gener- 
ally also appear at the low wavenumber end of the iner- 
tial subrange |2fi| | , and indications of this spectrum have 
been seen in convective dynamo simulations |27| . It turns 
out, however, that the value of C (obtained from a fit) is 
different for different values of Bq, casting doubt on the 
validity of the assumptions behind Eq. © . We therefore 
discard this simple explanation of the large scale mag- 
netic spectrum. Indeed in Fig. [7] we see that we get no 
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FIG. 8: Magnetic and kinetic power spectra for runs with 
S = 0.3 (runs B4, C4 and D4). 
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FIG. 10: Same as Fig. EJbut with B = 3.0. 
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FIG. 9: Fourier spectra of time evolution of the magnetic 
and kinetic fields at one point in the box for simulation with 
Bo = 0.06. The point of interest is chosen to be in the center 
of the box. The arrows represent the frequency of an Alfven 
wave with a wavelength of the box size traveling along the 
imposed field. We clearly see that Alfven waves are strongly 
present in the simulation. 



longer a fc _1 magnetic energy spectrum for large scales 
when the forcing is at k = 5; instead the infrared part of 
the spectrum has an increasing slope close to Eyi{k) ~ k 
for k < fcf. Some intermediate behavior is seen when 
fef = 2. ..3; see Ref. Q where no fc _1 behavior was found. 



VI. DIRECT EVIDENCE FOR ALFVEN WAVES 

Finally we look at the frequency power spectrum cal- 
culated from the time series of the magnetic field and 
velocity at one point in the simulation box; see Figs. 
and ^3 As expected, the larger the imposed magnetic 
field, the faster does the field oscillate. The peaks in the 
power spectra for Bq = 3.0 and Bq = 0.06 correspond to 



the frequency of the corresponding Alfven wave, 

co = v A ki, where v A = B /y/p, p (9) 

is the Alfven speed. (In our case we have /Lto = po = 1.) 
When Bq is comparable to or less than B cq , the peaks in 
the spectra are no longer well pronounced. 

For strong fields, however, the Alfven peaks are seen 
quite clearly. It is conceivable that these Alfven waves 
are stochastically excited by the turbulence. This might 
be similar to the stochastic driving of acoustic waves in 
the solar convection zone 1231 ■ 



VII. CONCLUSION 

The present studies have shown that a uniformly im- 
posed magnetic field has two important effects on the 
magnetic field that is induced at finite wavenumbers 
(k 0). First, the magnetic field is slightly enhanced 
at and around the forcing wavenumber (corresponding 
to the energy carrying scale). Second, the magnetic field 
is quenched with increasing Bq at all larger wavenumbers 
corresponding to the inertial and diffusive subranges. 

The enhancement and suppression at the two different 
wavenumber ranges is associated with a corresponding 
wavenumber dependence of the work term, Bq ■ (u x j). 
The suppression of the magnetic field in the inertial range 
is quite opposite to the behavior without imposed field 
when there is instead a significant enhancement of the 
magnetic energy spectrum over the kinetic energy spec- 
trum. We therefore refer to this effect as a suppression 
of the dynamo by the imposed field. 

The suppression of dynamo activity might be a conse- 
quence of the tendency toward two-dimensionalization of 
the turbulence by the large scale field |29|. Such an ef- 
fect is well-known for low-ReM hydromagnetic turbulence 
|30| , and it is a mathematical theorem that there can be 
no dynamo action in two dimensions |3lj . Of course, 
the turbulence does not really become two-dimensional, 
but instead the correlation length along the field becomes 
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large. This type of anisotrop y is a crucial ingredient of 
the Goldreich-Sridhar theory |18|. 

The Goldreich-Sridhar theory also predicts that Alfven 
waves should be present in the system. This has been 
confirmed by inspecting velocities and magnetic fields at 
a single point in the middle of the simulation box. These 
Alfven waves have the expected frequency uja = vpjt\. 
Furthermore, we do not find that there is equipartition 
between magnetic and kinetic energy spectra in the in- 
ertial range for large imposed field strengths. The ab- 
sence of equipartition may be a consequence of the iner- 
tial range being still too short (or absent). In runs where 
Bq = -B C q, on the other hand, there is clear evidence that 
kinetic and magnetic energy spectra fall on top of each 
other throughout the dissipation subrange. This is also 
in agreement with earlier results of Cho and collabora- 
tors [3 , who considered the case where the imposed field 
had equipartition strength. 

Whether or not models with imposed field can repro- 
duce the situation in small sub-domains of simulations 
with no overall imposed field is still unclear. At first 
glance the answer seems to be no, because none of the 
simulations with imposed field have ever been able to 
produce super-equipartition in the inertial range, as it 
is seen in the nonhelical simulations without imposed 
field ®. However, the reason for this may well lie in 
the still insufficient resolution of the simulations with no 



imposed field - even though they do already have a res- 
olution of 1024 3 meshpoints. It is indeed possible that, 
even though the kinetic and magnetic energy spectra are 
approximately parallel to each other over a certain range 
of wavenumbers and offset by a factor of about 2.5, they 
may actually converge at still larger wavenumbers. Pre- 
liminary indications of this have now been seen in sim- 
ulations using hyperviscosity and hyper-resistivity with 
no imposed field. However, a general difficulty with hy- 
perviscosity and hyper-resistivity is that certain aspects 
of the physics of such systems are significantly modified 
|32^ . It is therefore equally important to assess the fea- 
tures that are likely not to be altered by this manipula- 
tion. A detailed discussion of this will be the subject of 
a forthcoming paper. 
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